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Abstract 

Granular damping devices constitute an emerging technology for the attenuation 
of vibrations based on the dissipative nature of particle collisions. We show 
that the performance of such devices is independent of the material properties 
of the particles for working conditions where damping is optimal. Even the 
suppression of a dissipation mode (collisional or frictional) is unable to alter 
the response. We explain this phenomenon in terms of the inelastic collapse of 
granular materials. These findings provide a crucial standpoint for the design of 
such devices in order to achieve the desired low maintenance feature that makes 
particle dampers particularly suitable to harsh environments. 
Keywords: Particle dampers, Granular materials, Vibration attenuation 



1. Introduction 

Granular dampers (or particle dampers, PDs) are devices aimed at the at- 
tenuation of mechanical vibrations by exploiting the dissipative character of 
the interaction between macroscopic particles. A PD consists in a number of 
particles (grains) enclosed in a receptacle that is attached or embedded in a 
vibrating structure (see Fig. [T]). The motion of the grains inside the enclosure, 
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as the structure vibrates, is able to dissipate part of the energy through the 
non-conservative collisions, so reducing the vibration amplitude. This emerg- 
ing technology can replace the widely used viscous and viscoelastic dampers in 
particular applications where extreme temperatures (either low or high) are in- 
volved or where low maintenance is required. The leading sector in this regard 

nri n 

is the aerospace industry [!H3[. However, the automotive [J] and oil and gas 



industries are catching up in recent years. PDs are the descendant of the 



older impact dampers designed to operate by the use of a single body inside an 
enclosure [(j, 7|. PDs are now preferred over impact dampers due to the lower 
noise level they produce. The performance of a PD depends on a number of 
design characteristics such as the relative size and shape of the particles and 
the enclosure, the total weight m p of the particles, the number N of grains, 
the working vibration intensity and frequency, etc. These have been studied to 
some extent in the last two decades (8l4l0l|. However, less attention has been 
paid to the role that the grain-grain interaction plays in these systems. 

In this paper, we show that the response of a PD is independent of the 
particle-particle interaction to the extent that even friction or inelastic colli- 
sions can be suppressed without altering the vibration attenuation. This effect 
is explained in terms of the effective zero restitution of the granular system 
caused by the effective inelastic collapse 11 1. The inelastic collapse in dense 



granular materials refers to the effect by which the system can dissipate its entire 
kinetic energy in a short finite time even if collisions have a very high restitution 
coefficient. In dense systems, the number of collisions grows so rapidly that even 
a minute dissipation in each collision suffices to make the system as a whole fully 
dissipative. We show that this interpretation allows us to set the limits to the 
universal response. Crucial implications for the design and maintenance of PDs 
are discussed. 
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Figure 1: (Color online) Schematic view of a particle damper. (M) the primary mass of the 
structure, (m p ) the total mass of the particles in the enclosure, (K) the spring constant, (C) 
the viscous structural damping, (B) the vibrating base where the displacement is imposed. 

2. DEM Simulations 

Wc have carried out simulations of a PD by means of a Discrete Element 
Method (DEM). A number N (with 5 < N < 250) of spherical particles of 
mass to = 9.08 x 10~ 4 kg are deposited in a prismatic enclosure (see Fig. [IJ 
of mass M = 2.37 kg, base L X L (L = 0.03675 m) and height L z (with 0.057 
m< L z < 0.282 m). The enclosure is attached to a vibrating base by means 
of a Hookean spring (spring constant K — 21500 Nm _1 ) and a viscous damper 
of low dissipation constant (C = 7.6 Nsm -1 ). The particles interact through a 
Hertz-Kuwabara-Kono [l2j normal contact force (F n — —k n a 3 / 2 —j n v n y/a) plus 
a tangential force (Ft — — mm (|7tUt\/a| ,\jUdF n \) sgn (v t )) that implements the 
frictional property of the grain surfaces [13| |. Here, v n and Vt are the relative 



3 



normal and tangential velocities and a = rij — d the virtual overlap between the 
interacting particles i and j. The diameter of the particles is d = 0.006m. The 
parameters that set the grain-grain interactions are: k n = 4.02 x 10 9 Nm~ 3 / 2 
(which corresponds to the Young modulus E — 2.03 x 10 11 Nm" 2 and Poi- 
son ratio v = 0.28 for steel), < j n < 1 x 10 4 kgm _1 / 2 s _1 , < 7 f < 1 x 10 4 
kgm _1 / 2 s~ 1 and /^d = 0.3. The grain- wall interaction is taken equal to the grain- 
grain interaction. The vibrating base and the enclosure are constrained to move 
only in the vertical z-direction, whereas the particles can move in the three- 
dimensional volume of the receptacle subjected to the action of gravity. The mo- 
tion Zbase(t) of the base is set to a harmonic function [zb a se(t) = 0.0045 cos(wt) 
m] whose frequency to — 2nf is varied in the range 0.5Hz < / < 30Hz. We 
monitor the amplitude of the oscillation z max of the enclosure in response to 
the base vibration, the motion of the grains inside and the energy dissipated in 
each cycle. All analysis is done over the stationary state of the system for each 
given / after any transient has died out. 

Figure [2f a) shows the frequency response function (FRF) — i.e., z max as a 
function of / — of the system with N = 250 for a given choice of the gran- 
ular interaction and three different height L z of the enclosure in comparison 
with the response obtained when the enclosure is empty. As discussed in pre- 
vious studies 0, Q, Q| ) there exists an optimum enclosure height (in this case 
L z — 0.12255m) for which the maximum attenuation is achieved. The resonant 
frequency is shifted due to the added mass (the granular mass), but in a non- 
trivial way, with overshoot and undershoot effective masses [it] ] . The motion of 
the granular bed inside the enclosure for the optimum enclosure height can be 
seen in Fig. EJb), (c) and (d) for different values of /. The grains behave as a 
more or less dense lumped mass for a wide range of frequencies (0 < / < 18Hz). 
However, for high frequencies and tall enclosures, the grains enter a gas-like 
state 171. 
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3. Effect of particle particle interactions 

In order to asses the effect of the particle-particle interactions, we plot in 
Fig. [3ja) the FRF for the optimum enclosure height for different values of 
7„ and j t - The effective normal and tangential restitution coefficient is an 
exponentially decaying function of these parameters and depend on the relative 
velocity at impact [13j . As we can see, different interaction parameters yield the 
same FRF, suggesting a universal response. Notice that even eliminating the 
frictional character of the particles or, alternatively, eliminating the dissipative 
nature of the normal interaction is not sufficient to induce a change in the FRF 
over a wide range of frequencies. As it is to be expected, eliminating both, the 
normal and tangential dissipative part of the interaction converts the system 
into a conservative molecular-like system which yields no attenuation of the 
response [see diamonds in Fig. Eta)]. In this case only a shift of the resonant 
frequency is observed due to the added mass m p [see brown dashed line in Fig. 
Elja) for the analytical solution]. We have also considered different material 
properties such as Young modulus and dynamic friction coefficient, but have 
observed no change in the FRF. 

The energy dissipated in each oscillation cycle is shown in Fig. E{b) as 
a function of the excitation frequency. Notice that the frictional dissipation 
and the collisional dissipation are higher around the resonant frequency. The 
proportion of energy dissipated through one mode or the other (collisional or 
frictional) depends on the actual interaction parameters. However, the total 
energy dissipated is independent of the friction and restitution characteristics 
of the particles. 

To further explore the extent of the apparent universal response of the PD, 
we have considered a wide range of the dissipative interaction parameters 7„ 
and 7 t . As it can be seen in Figs. Etc), where the response of the system with 
N = 250 at the resonant frequency is considered, increasing 7„ at constant 
7 t leads to an increase in the energy dissipated though the collisional mode. 
The converse is true if the frictional parameter j t is increased [see Fig. EJd)]. 
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However, the total dissipated energy remains constant even if j n or 74 drops 
to zero. Hence, the system is able to dissipate the same amount of energy 
irrespective of the dissipation modes available to the particles. 



4. Origin of the universal response 

We speculate that this universal response is found whenever a large number 
of particles is used and the motion of the granular bed is set into a more or less 
dense lumped mass, as oppose to a gas-like state. For a dense granular layer, 
the number of collisions per unit time as the bed collides with the boundaries 
increases dramatically due to an effective inelastic collapse . Although 

inelastic collapse refers to a mathematical divergence of the number of collisions 
per unit time when instantaneous interactions are considered, real systems do 
also exhibit a remarkable increases in the collision rate . The behavior of 

such granular systems has been recently proven to be well modeled by a single 
mass m with an effective zero restitution coefficient and no frictional properties 
[2c] ]. Also, recent simulations with no frictional components where able to fit 
experimental data on PDs in microgravity [21| . This fact — that friction can 
be neglected — is evidenced in Fig. Etd). However, according to Fig. Etc), 
restitution can also be set to unity and any non- vanishing friction will suffice to 
render the universal FRF shown in Fig. Eta). This indicates that a much wider 
set of interactions can lead to the effective inelastic collapse than previously 
shown. 

We have solved a simple one-dimensional model where a single particle of 
mass nip and zero restitution coefficient moves between the floor an ceiling. 
The model is adapted from [8(. The results are shown in Figs. Eta) andEtb) 
with black solid lines. It is clear that such simple model provides the essence 
to describe the PD. Hence, design can be based in this simple model without 
worrying about a careful selection of the particle properties. Previous workers 
have used slightly more complicated models where the restitution coefficient was 

Our results indicate that this complication 



used as a fitting parameter [S 
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may be unnecessary. 

5. Limits to the universal response 

We now turn into surveying the limits of this universal response. Inelastic 
collapse is known to take place only if a large number of grains at high densities 
are considered. Therefore, we can set the system into a regime where this 
universal response does not apply by reducing TV or by promoting a diluted 
granular state in the enclosure. 

In Fig. 21 we present the variation of the amplitude of vibration z r at the 
resonant frequency as a function of the normal [Fig. Ufa)] and tangential [Fig. 
0Jb)] dissipative parameters for different N. For these simulations, we have 
changed the diameter of the particles in order to keep the total mass m p of the 
grains constant as TV is changed. For N > 100, a constant {universal) response 
is recovered, whereas smaller systems present a better attenuation as any of 
the two dissipative properties (p/ n or j t ) is increased, in accord with intuition. 
These results confirm the speculation that the universal response only applies 
if a relatively large number of particles is involved. Of course, the values of N 
at which this universal response is reached will depend on the horizontal size of 
the enclosure. We estimate from our simulations that whenever the enclosure is 
filled with three or more layers of particles the system response near resonance 
becomes independent of the particle-particle interaction. 

A way to induce a gas-like behavior of the granular sample in the enclosure 
— so as to create a dilute regime where an effective inelastic collapse is not 
expected — is to increase the height of the cavity. In Fig. EJa) and (b) we 
plot the energy dissipated per cycle as a function of 7„ and 7 t for an enclosure 
with L z = 0.282m. For this large L z , the granular sample expands significantly 
and does not move as a lumped mass, so reducing dramatically the number 
of collisions per unit time. The final results is an effective dissipation that 
depends on the particle-particle dissipative interaction. Interestingly, in this 
regime, increasing j n or j t leads to a decrease of the total energy dissipated. 
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Since PDs are designed to optimize attenuation in a number of applications, the 
size of the enclosure generally promotes the dense lumped mass motion of the 
grains inside [1]. Therefore, in many working conditions of interest, the system 
will be in a regime where the universal FRF is obtained. 



6. Conclusions 

We have shown that a basic phenomenon (inelastic collapse) leads to a uni- 
versal response of a PD — in the sense that the particle-particle interaction 
becomes irrelevant. This allowed us to determine the limits of this universality: 
relatively large numbers of particles must be used and the system has to be set 
in a state of dense lumped mass. 

It is worth mentioning that the particular contact force model used for the 
simulations is of little relevance to the results presented here. We have shown 
that a much simplistic model (the zero-restitution-singlc-particle model) also 
shows the same response. Indeed, the proposed universality implies that details 
of the interactions are not relevant. Moreover, despite our studies being done 
on a system vibrating in the direction of gravity, we expect results will apply to 
horizontally vibrating PDs. 

This universal response is consistent with some observations in experiments 
and simulations where a few values of the material properties where tested 
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261 ] and can explain the unexpected agreement between simplified models and 
complex experiments [21| . It is worth mentioning that powders, as opposed to 
granulars, may not follow this universal response even at resonance. Powders are 
fine graded particles and the effects of the hydrodynamics of the surrounding air 
affects the motion of the particles to a large extent. Fine powders will expand 
due to air-particle interactions and the inelastic collapse will be unlikely. Some 



preliminary experiments with PDs using fine powders seem to confirm this 26]. 

The suitability of PDs to work in harsh environments can be understood as 
a consequence of this phenomenon. Extreme temperatures and pressures may 
induce mild changes in frictional properties, but these will not alter the PD 
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response. More importantly, degradation of the particles during operation due 
to wear, deformation and fragmentation are not likely to compromise the PD 
performance. Changes in friction or restitution are unimportant. Although we 
have not studied particles of different shapes, we speculate that fragments of 
particles may be as effective as the original particles as long as they are not fine 
graded. Moreover, fragmentation can only increase N, which should not take 
the system out of the universal FRF. This is the underlying phenomenon that 
explains the characteristic low maintenance required for this devices. Notice 
however that very high temperatures may weld particles together inducing an 
effective reduction of N which can reduce vibration attenuation. 

Design of PDs can be greatly simplified by choosing to work with large N and 
using a simple model such as the zero restitution single mass used here. Under 
these conditions, the selection of the particle material properties is unimportant 
for the PD performance and one can focus, for example, on cost effectiveness. 
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Figure Captions 

Figure 1: (Color online) Schematic view of a particle damper. (M) the pri- 
mary mass of the structure, (m p ) the total mass of the particles in the enclosure, 
(K) the spring constant, (C) the viscous structural damping, (B) the vibrating 
base where the displacement is imposed. 

Figure 2: (Color online) (a) Frequency response function of the system with 
N = 250, 7„ = 3660.0 kgrn" 1 /^- 1 and 7t = 10980.0 kgrn" 1 /^- 1 . We plot the 
maximum amplitude of the oscillation z max of the primary mass as a function 
of the excitation frequency / imposed to the base for different heights L z of 
the enclosure (see legend) and compare with the theoretical FRF for an empty 
enclosure (black solid line). Lines joining symbols are only to guide the eye. 
In panels (b), (c) and (d), we plot the trajectory of the enclosure (black solid 
lines define the floor and ceiling) and the motion of the granular sample inside 
(colored area defined by the position of the uppermost and lowermost particle 
at each time; the dotted red line is the position of the center of mass) for the 
optimum height L z = 0.12255m: (b) / = 5.5Hz, (c) 14.5Hz and (d) / = 21.0Hz. 

Figure 3: (Color online) (a) The FRF for a PD of optimum height (L z = 
0.12255m) for N = 250 particles with different interaction parameters 7 „ and 7t 
(see legend for the values used in units of kgm _1 / 2 s _1 ). The dashed lines corre- 
spond to the analytical solution for: an empty enclosure (orange), and an empty 
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enclosure with primary mass corrected to M + m p (brown). The black solid line 
corresponds to the zero restitution single particle model, (b) The energy dis- 
sipated per cycle as a function of / for the frictional (squares) and collisional 
(triangles) modes for three different particle-particle interactions (see legend). 
The total energy dissipated (circles) is independent of 7„ and 74. The black solid 
line corresponds to the zero restitution single particle model, (c) Energy dissi- 
pated as a function of the collisional dissipation 7„ at the resonant frequency for 
7i = 10980.0 kgm _1 / 2 s _1 . (d) Energy dissipated as a function of the frictional 
dissipation j t at the resonant frequency for 7„ = 3660.0 kgm _1 / 2 s _1 . 

Figure 4: (Color online) (a) The vibration amplitude z r at resonance as 
a function of the collisional dissipation j n for j t — 10980.0 kgm _1 / 2 s _1 for 
different number N of grains in the enclosure (see legend). The particle size is 
chosen to yield a total particle mass m p = 0.227 kg. (b) Same as (a) but the 
dependence on j t is considered for 7„ = 3660.0 kgm _1 / 2 s _1 . 

Figure 5: (Color online) (a) The energy dissipated per cycle at resonance as 
a function of the collisional dissipation -f n for a large enclosure (L z = 0.282m) 
with N = 250 particles and 7* = 10980.0 kgm" 1 /^- 1 . (b) Same as (a) but the 
effect of the frictional dissipation 74 is considered for j n = 3660.0 kgm _1 / 2 s _1 . 



13 



0.16 
0.12 

' 0.08 

j 

0.04 
0.00 



Lz = 0.057 m 


(a) 


b Lz = 0.12255 m n 




" e Lz = 0.282 m « 




- Empty enclosure 









-0.06 




-0.06 



5.0 7.5 10.0 12.5 15.0 17.5 20.0 
/ [Hz] 



11.30 11.41 




-0.06 



11.30 11.41 



11.51 11.62 
t [S] 



11.72 



11.51 11.62 

t rsi 



11.72 



(d) 



/=21.0 Hz 



11.30 11.41 



11.51 11.62 
t [S] 



11.72 



Figure 2: (Color online) (a) Frequency response function of the system with N = 250, 7n = 
3660.0 kgm - 1 / 2 s _1 and 7t = 10980.0 kgm~ 1 / 2 s~ 1 . We plot the maximum amplitude of the 
oscillation z max of the primary mass as a function of the excitation frequency / imposed to 
the base for different heights L z of the enclosure (see legend) and compare with the theoretical 
FRF for an empty enclosure (black solid line). Lines joining symbols are only to guide the 
eye. In panels (b), (c) and (d), we plot the trajectory of the enclosure (black solid lines define 
the floor and ceiling) and the motion of the granular sample inside (colored area defined by 
the position of the uppermost and lowermost particle at each time; the dotted red line is the 
position of the center of mass) for the optimum height L z = 0.12255m: (b) / = 5.5Hz, (c) 
14.5Hz and (d) / = 21.0Hz. 
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Figure 3: (Color online) (a) The FRF for a PD of optimum height (L z = 0.12255m) for 
N = 250 particles with different interaction parameters 7„ and -ft (see legend for the values 
used in units of kgm _1 / 2 s — 1 ). The dashed lines correspond to the analytical solution for: an 
empty enclosure (orange), and an empty enclosure with primary mass corrected to M + m p 
(brown). The black solid line corresponds to the zero restitution single particle model, (b) 
The energy dissipated per cycle as a function of / for the frictional (squares) and collisional 
(triangles) modes for three different particle-particle interactions (see legend). The total en- 
ergy dissipated (circles) is independent of 7„ and 74. The black solid line corresponds to the 
zero restitution single particle model, (c) Energy dissipated as a function of the collisional 
dissipation j n at the resonant frequency for 74 = 10980.0 kgm - 1 / 2 s — 1 . (d) Energy dissi- 
pated as a function of the frictional dissipation 74 at the resonant frequency for 7„ = 3660.0 
kgm - 1//2 s — 1 . 
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Figure 4: (Color online) (a) The vibration amplitude z r at resonance as a function of the 
collisional dissipation 7„ for yt = 10980.0 kgm~ 1//2 s — 1 for different number TV of grains in the 
enclosure (see legend). The particle size is chosen to yield a total particle mass m p = 0.227 
kg. (b) Same as (a) but the dependence on jt is considered for j n = 3660.0 kgm _1 / 2 s — 1 . 




Figure 5: (Color online) (a) The energy dissipated per cycle at resonance as a function of 
the collisional dissipation 7 n for a large enclosure [h z = 0.282m) with N = 250 particles and 
"/t = 10980.0 kgm - -^s - 1 . (b) Same as (a) but the effect of the frictional dissipation ft is 
considered for 7„ = 3660.0 kgm _1 / 2 s _1 . 
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